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ABSTRACT 

Recent calculations of the hydromagnetic deformation of a stratified, non-barotropic 
neutron star are generalized to describe objects with superconducting interiors, whose 
magnetic permeability fi is much smaller than the vacuum value ^o- It is found that 
the star remains oblate if the poloidal magnetic field energy is > 40% of total mag- 
netic field energy, that the toroidal field is confined to a torus which shrinks as fi 
decreases, and that the deformation is much larger (by a factor ~ Mo/m) than in a 
non-superconducting object. The results are applied to the latest direct and indirect 
upper limits on gravitational-wave emission from Laser Interferometer Gravitational 
Wave Observatory (LIGO) and radio pulse timing (spin-down) observations of 81 mil- 
lisecond pulsars, to show how one can use these observations to infer the internal field 
strength. It is found that the indirect spin-down limits already imply astrophysically 
interesting constraints on the poloidal-toroidal field ratio and diamagnetic shielding 
factor (by which accretion reduces the observable external magnetic field, e.g. by 
burial). These constraints will improve following gravitational- wave detections, with 
implications for accretion-driven magnetic field evolution in recycled pulsars and the 
hydromagnetic stability of these objects' interiors. 
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1 INTRODUCTION 



The external magnetic field of a neutron star is (relatively) easily inferred from its spin-down rate, but its internal magnetic 
field is not directly observable. The main clue suggesting the existence of strong in ternal neutron star fields comes from 
the 1998 August 2 7 giant flare from the soft gamma-ray repeater (SGR) 1900+14 jFeroci et al.lll999l : lliurlev et al J 1 19991 : 



Mazets et al 



19991 ) . The giant flare, which released ~ 10 3 ^ J of energy as X-ray s, was accompani ed by a 2. 3- fold increase in 



the spin-down rate ( Mazets et al.|[l99^ : Iwbods et al. 19991 : Thompson et al. 2000). T° explain this. Ilokal (2001) proposed that 
the flare and the enhanced spin down were cau sed by a global reconf iguration of the internal ma g netic field of ~ 10 13 T, well 
above the external dipole field of 6.4 x 10 10 ifl lCorsi fc Owenl ifcoill ) generalised the Ilokal (120011 ) calculation (by allowing the 
toroidal field strength to change, as well as the moment of inertia) and c oncluded th at an internal field strength of ~ 10 12 
T gave rise to the 1998 August 27 event, lower than the first estimate by Ilokal (|200lr ). but still significantly higher than the 
observed external field. 



Stellar ellipticity can also b e used to constrain the strength of a star's internal field (I Cutler 2002 



Abbott et al 



201 



DalPOsso et al 



2009 



Pitkin l201ll). It is well known that a strong magnetic field can deform a star Jchandrasekhar fc Fermi 



1953l : lFerrarolll954l : lGoosenslll972l : lKatj|l989l : bavne fc Melatos|l2004lHaskell et ap2008l;lMastrano et alJboill ). The ellipticity 



e is roughly proportional to the magnetic energy ( Cii7leiTl2'u02nHask6^ Neutron stars, with 
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their intense magnetic fields, possess significant elliptic i ties, making them good candidates for gravitational wave sources 
( Bonazzola fc Gourgoulhonlll996l ; lMelatos fc Pavnell2005l ; Istella et all2005l ; lHaskell et alJliooi i lDalPOsso et a,l.|[20Qg| v Recent 
data from the fifth Laser Interfe rometer Gravitatio n al W ave Observatory (LIGO) Science Run set an upper limit of e < 



1.4 x 10 on the Crab Pulsar (| Abbott et al.ll200Sl , |20ich . translating into an internal magnetic field of < 10 12 T under 



standard assumptions. LIGO non-detections of the central compact object (CCO) in the supernova remnant Cassiopeia A 
(Cas A) have constrained its ellipticity as well. The Cas A CCO has not been detected electr omag netically. mak ing it impossible 
to infer its external magnetic field from the spin-down rate. However, Wette et al. 1 20081 ) and Wette! ( 2010h constrained its 
ellipticity as a function of gravitational wave frequency (e.g., e < 3.6 x 10" 4 fo r 100 Hz, e < 0.6 x 10" 4 for 200 Hz, and 
4 for 300 Hz), implying an internal magnetic field < 10 14 T. Lastly. IChung et al.l ( 201 If ) showed that it will be 



of the putative 24-year-old 



e < 0.38 x 10" 

possible to use future data from LIGO to set a lower limit of 10 7 T on the magnetic field (e < 10" 
neutron star in the supernova remnant SNR 1987A. 

Gravitational waves are generated by a rotating star when it is not spherically symmetric and when its 'wobble angle' 8, 
the angle between its total angular momentum vector and symmetry axis, is nonzero. In general, therefore, the star precesses 
as it radiates. The magnetic field analysed in this paper is axisymmetric, which deforms the star into an ellipsoid. The 
most general expression for the grav itational wave signal of a precessing triaxial ellipsoid [given by, e.g., Eqs. (19)-(26) of 
Uaranowski. Krolak. fc Schutzll200ll )] depends on 9 and the angle i between the angular momentum vector and line of sight 
to the Earth, as well as e. The signal is strongest when 9 = 7r/2 (i.e., the rotation axis is perpendicular to the symmetry axis, 
which is in turn parallel to the magnetic axis) and l — (i.e., the rotation axis is directed towards Earth). In this paper, to 
simplify matters, we henceforth assume implicitly that 9 = tv/2 and i = 0, so that gravitational wave emission and signal 
detecti on are assumed t o be o ptimal and there is no precession, in order to focus on e without geometric complications. 

( 20111 ). we constructed hydromagnetic equilibria for stratified, non-barotropic stars. The commonly 



In Mastrano et al 



adopted bar otropic assumption, while simplifying calculations, severely restricts the form of the field that ca n be 'fitted' into 



the st ar, e.g., Haske ll et al.l (120 08) foun d that the field must vanish at the surface, contrary to observations, and lLander fc Jones 



( 2009) and lCiolfi. Ferrari, fc Gualtieril (|20ld ) found that only configurations dominated by the poloidal com ponent (poloidal 
energy > 90% of total) are allowed, contrary to the numerical simulations o f lBraithwaite fc Nordlundl 12006! ). By abandoning 
the barotropic assumption, we are able to construct a simple, self-consistent hydromagnetic equilibrium with an internal 
field that can be matched to an external dipole. In our configuration, because we do not require the pressure to be a 
function of density alone, we are less restricted in the choices of poloidal and toroidal components; they need not be of any 
particular relative strengths and are independentl y adjustable. Inc i denta lly, t his means that the magnetic field configuration 
is independent of the equation of state chosen, cf. Lander fc Jones! 1 20091 ) and Ciolfi. Ferrari, fc Gualtieri |2010l ). 

In this paper, we present one possible astrophysical application of the aforementioned result, namely to constrain the 
internal magnetic fields of neutron stars in conjunction with gravitational wave observations. As we can match the internal 
field to an external dipole, we are able to relate one set of observations (external magnetic field strength, from spin period 
and spin-down rate) with another (gravitat ional wave upper limit s), at least in principle. In Sec. 2, we briefly describe our 
field structure and summarize the results of iMastrano et alj (|201ll ). Then, in Sec. 3, we generalize the earlier work to include 
a superconducting interior. We compare the ellipticity calculated using our superconducting model to the gravitational-wave 
upper limits of 81 known millisecond pulsars, including also the important effect of accretion- induced diamagnetic shielding. 
Lastly, in Sec. 4, we summarize our results and discuss the possibility that the internal fields of millisecond pulsars are stronger 
than currently thought. 



2 FIELD STRENGTH VERSUS ELLIPTICITY FOR NON-BAROTROPIC STARS 



Mastrano et al.l (|201ll ) considered a general class of poloidal-toroi dal magnetic field configurations, which are broadly repre- 
sentative of the fi eld structures observed in numerical simulations ( Braithwaite fc Nordhmd 20061 ; Braithwaite fc Spruit 20061 ; 
Braithwaite1l2009l ) and exhibit the following properties: 



• the field is axially symmetric around the z-axis; 

• the poloidal part is continuous with a dipole field outside the star (so there are no surface currents) and vanishes at 
the circle (r, 9) — (0.78-R, 7r/2) (where r and 9 are the radial and polar coordinates and R is the stellar radius; this 
locus is called the neutral circle); 

• the toroidal component is confined to the region of closed poloidal field lines around the neutral circle; 

• the current density remains finite and continuous everywhere in the star. 

We write the magnetic field in the form pioneered by Chandrasekhar 1 19561) , 



n = B [n P Va(r,9)xV<l> + r]tP(a)V<f>], (1) 
where r\ v and rjt are dimensionless parameters which define the relative strengths of the poloidal and toroidal components 
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respectively. The function /3(a) takes the form /3(a) = (a — l) 2 for a ^ 1 and /3(a) = elsewhere, confining the toroidal field 
to the region where a exceeds unity, the value taken by a at (r, 8) = (1, 7r/2); the current density goes continuously to zero at 
this boundary. The flux function a(r, 8) is taken to be f(r) sin 2 8. Note that this particular form of a(r, 8) is only applicable 



when we try to match our field to an external dipole; other multipoles match different 
of a is given by 



-dependences. The radial dependence 



m = f 



6r 4 3r' 



(2) 



The function f(r) is postulated to be of this form to ensure that the field described by Eqs. (l)-(2) is continuous with a 
dipole field outside the s tar, that ther e are n o surface currents, and that the current density is finite at the origin [for a more 
thorough derivation, see lAkgiin et all |201lh ]. 



A schematic diagram of the field is shown in Fig. [11 I It represents just one possibility amongst many, chosen for simplicity 



and mathematical convenience; possible conf igurations involving higher multipoles are not ruled out by observations ( Aronsl 
19931 ; iThompson. Lvutikov. fc Kulkarnill2002h . 

We now calculate the small changes Sp and Sp to the pressure and den s ity of a star in hydrostatic equilibrium caused by 
this field, which satisfy the force balance equation ( Easson fe Pethick 1977 ; Akgiin fc Wasserman 20081 ) 



(V x H) x B = VSp + SpV<f>. (3) 

In Eq. (3), H = B//i is the magnetic intensity, $ is the gravitational potential (the Cowling approximation has been taken, 
with 5$ = 0), and p is the magnetic permeability. This allows for a superconducting interior (where p is smaller than the 
vacuum permeability /Uo), as analysed in Sec. 3. The assumption that the changes to density and pressure are small enough 
that they can be treated as perturbations on the steady state is justified a posteriori. Ellipticity is then calculated from 
the perturbed density. Again, the non-barotropic assumption is essential. If barotropy is assumed instead, pressure must 
be a function of density only, to all orders; a restriction is then imposed on the magnetic field confi guration, because th e 
poloidal and tor o idal c ompo nents must be related in su c h a w ay that Sp is purely a function of 8p [e.g.. Irlaskell et al ( 20081 ) . 
lLander fc Jonesl l|2009l ). and lCiolfi. Ferrari, fc Gualtieril |20ld )]. 

For the initial unperturbed hydrostatic equilibrium, we adopt a parabolic density profile 



p = p c (l-r 2 ), (4) 

where p c = 15M*/(87ri?2) is the density at the core, M* is the mass of the star, and i?* is the radius of the star. While 
this is one part i cular simple choice of density profile, chosen to render the calculations tractable, we showed in Sec. 3.3 of 
Mastrano et al.l l|201ll ) that the resulting ellipticity is within 5% of that obtained using the more common n = 1 polytrope. 

We define the parameter A as the ratio of poloidal field energy to total field energy. The energies are obtained by 
integrating the squa res of the poloidal and total magnetic intensities over the star; note that this definition is slightly different 
from that given by iMastrano et al.l (201l|), who integrated H 2 over all space. We have A = 1 for a purely poloidal field 
configuration and A = for a purely toroidal configuration. In terms of A, for a non-superconducting star with p = po, the 
ellipticity e takes the form 



M, 



( R* 



0.351 



where B B is the surface magnetic field strength at the equator, and Mq is the Solar mass. As expected, one finds e oc B B . 
The mass quadrupole moment vanishes for A = 0.351, that is, when the poloidal f ield energy is 35.1% of the t otal m agnetic 
field energy. Note that both our model and the generalised 'twisted torus' model of ICiolfi. Ferrari, fc Gualtieri |2010h predict 



e ~ 4 x 10 6 for a canonical 'magnetar-like' neutron star with a purely poloidal field of strength B B = 5 x 10 10 T. 



3 MILLISECOND PULSARS 



Armed with Eq. (5), we can use spin-down measurements of B s and gravitational- wave upper limits on e to constrain A and 
hence the internal magnetic field for various classes of object. 

In Fig. [2j we plot the most up-to-date LIGO (triangles) and spin-down (dots) gravitationa l-wave upper limits on the 
absolute value of e for 81 known millisecond radio pulsars as a function of B B ( Abbott et al. 2Q10l) . As expected, the indirect 
spin-down limits are uniformly tighter than the direct LIGO limits for now, although this will change in the future; already, 



2 The stability of this configuration is examined bv lAkgun et al. 1 1 120111) . 
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Figure 1. A cross-sectional diagram of the magnetic field structure analysed in this paper. The surface of the star is represented by the 
dashed circle. The toroidal magnetic component is confined to the shaded region and fills a torus around the z-axis. The poloidal field 
lines vanish at the neutral circle, i.e. at the dots in the shaded regions. 



LIGO has be aten the spin-down limits for a handful of non-millisecond pulsars like the Crab ( Abbott et al. 2008) and PSR 
J0537-6910 (| Abbott et allboioh . However, the LIGO limits are included in Fig. [2] to show, in the discussion below, what 
constraints can be extracted from gravitational-wave detections when a spin-down measurement is not available, e.g. for a 
CCO like Cas A. Fig. [2] also displays a selection of theoretical curves e(-B s ), seen as diagonal lines in the figure. The curves 
encompass a range of evolutionary scenarios discussed in turn in the subsections below. 

Let us begin by considering the lower black dashed curve, which corresponds to Eq. (5) with fi = fio and A = 10 -3 . 
In this scenario, the star is not superconducting, the external field is not reduced by accretion in any way, and the internal 
toroidal field is relatively strong, with the poloidal component contributing only 10 -3 to the total magnetic energy. Even so, 
e lies ~ 6 orders of magnitude below the strictest of the spin-down limits. 

Next, consider the upper black dashed curve, which corresponds to Eq. (5) with fi = fio and A = 10 -10 . Now the theoretical 
curve sits comfortably in the region spanned by the spin-down data points. However, this curve represents a physically extreme 
case, where the toroidal field strength is ~ 1% of the virial field [~ 10 14 (M*/1.4M o )(i?,,/10 4 m)~ 2 T jLattimer fc Prakash 
2007)] and the toroid al field energy is ~ 10 10 times the poloidal field energy. The value of A is far below the stable lower limit 
of 10 -3 calculated by Braithwaite ( 20091 ) . suggesting that the upper black dashed curve cannot be realised in practice. 

Two important extra effects act in concert to bring the theoretical curves close to the data without appealing to extreme 
situations like the one in the previous paragraph. First, all objects in Fig. [2] are intentionally selected to be recycled. It is 
likely that the actual internal magnetic field is much stronger than the measu red external dipole field in a recycled pulsar, 
because surface currents are diamagnetically shielded ( Bisnovatvi-Koean 1974), buried by polar accretion ( Pavne fc Melatoa 
20041 : IMelatos fc Pavne] l2005h. or resistively dissipated J Romanil 119901). Second, it is also likely that the core of t he star is 
a superconductor |joneslll975l ; lEasson fc Pethidjfl977l : ICutlerll2002l : IWassermanlboOoj ; lAkgiin fc Wassermanl [20081 ). We now 
discuss these scenarios in turn. 



3.1 Core superconductivity 



If the stellar core \p c > 2.8 x 10 17 kg m 3 iLevenfish. Shibanov. fc Yakovlev 199S )] is made of a type II superconductor 
llBaym. Pethick. fc Pinesll969l ; lBavm fc Pethickll975l ; lEasson fc Pethickll977l:lElgarOv et al.ll99r3:|jonesl20oilBaldo fc Schulzj 



ivm, Pethick. & Fmesll9rJ9l : IBavm & Fethickll975l ; IEasson & Feth ick 1977; Elgam 
2007l h the magnetic permeability fi ~ 10 3 /jq drops significantly faasson fc Pethicklll977l : lAkgiin fc Wassermanll2008l ) and Sp 



and 5p are magnified by the factor ^xn /^t. lAkgiin fc Wassermanl |2008l ) constructed a model of such a star, where the supercon- 
ducting core is surrounded by a Newtonian fluid, but their model assumes that the star is barotropic and the internal field is 
purely toroidal. We can generalize their result to our particular non-barotropic configuration. For simplicity, we assume that the 
entir e star is superconducting in this first attempt, instead of just a core region up to 0.6i?, jYakovlev. Levenfish. fc Shibanovl 
19991 ). 
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Figure 2. Direct LIGO upper limits (triangles) and indirect radio timing spin-down limits (dot s) on the absolute ellipticities |e| of 
81 known millisecond pulsars as a function of their measured surface magnetic field strengths B s \ Abbott et al.ll201Cf) . Also shown are 
theoretical |e(Z? s )| curves for the following structural and evolutionary scenarios: non-superconducting (/i = fig) (the dashed black curves), 
and superconducting (fi = 10~ 3 /io) with diamagnetic shielding (red and blue curves). The non-superconducting curves are for A = 10~ 3 
(realistic; lower dashed black curve) and A = 10~ 10 (unrealistic; upper dashed black curve). The red and blue bands correspond to 
diamagnetic shielding factors £ = B s observed/^ s actual = 10 — 4 and 10 -2 respectively. Specifically we have A = 10 -3 , § = 10 -4 (upper 
solid red curve), A = 1, § = 10 -4 (lower solid red curve), A = 10 — 3 , £ = 10 — 2 (upper solid blue curve), and A = 1, § = 10~ 2 (lower solid 
blue curve). Note that all the A = 1 curves correspond to oblate stars and all the A = 10 curves correspond to prolate stars. 

Firstly, to ensure that the stellar field still possesses the propertied outlined in Sec. 2, the function /(r) in Eq. (2) must 
be modified into 

f(r) = -[(26 + 9^)r 2 - (28 + U^ r )r 4 + (10 + 5^-)r 6 ], (6) 
o 

with n r = fi/fio [note that this reduces to Eq. (2) for fi r = 1]. The magnetic field is modified in two ways: (1) the field lines 
become more radial just inside the surface; and (2) the region which contains the toroidal field squeezes closer to the surface 
and shrinks. In Fig. [3] we sketch the field lines for fi r = 10 _1 [Fig. EJa)] and 10 -3 [Fig. [Sfc)]. Recall that the toroidal field is 
confined to the region where f(r) sin 2 8^1 (this is also the region described by the outermost poloidal field line which closes 
inside the star). As evident from Eq. (6), this region shrinks as \x r decreases (but does not vanish for ji r = 0, the case of 
perfect superconductivity, tending instead to ~ 7.5 x 10~ 3 of the total volume in the limit fi r — > 0). We magnify this region 
in Fig. E^b) (for /i r = 1CP 1 ) and (d) (for fi r = 10 -3 ). The volume of the torus is 0.16 (0.04) of the original /i r = 1 case for 
fir = 10" 1 (fj, r = 10" 3 ). 

The magnetic intensity H = B//x in Eq. (3) is magnified by n/no in a superconductor. The deformation caused by this 



field is then calculated by the method given by iMastrano et al.l l|201ll ). We are unable to derive a simple analytic formula 



3 Specifically, we must ensure that the normal component of the internal field and the tangential component of the magnetic intensity 
are continuous with the external dipole field at r = 1. 
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Table 1. Dimensionless constants C12 in Eq. 0, for some selected values of /i r = fi/no- 



Mr 


Cl 




C2 


10° 


5.62 x 10' 


6 


0.35 


5 x ict 1 


9.97 x 10- 


-6 


0.42 


10- 1 


3.92 x 1Q- 


-5 


0.40 


HT 2 


3.84 x 10" 


-4 


0.42 


10- 3 


3.83 x 10" 


-3 


0.43 


10- 4 


3.84 x 10" 


-2 


0.43 



describing the dependence of e on fi r = ft/fta and A simultaneously, but we find that the general form of e(A) for a given fj, r 
is still similar to Eq. (5), namely 

c = Ci ( 5 x101°t) (l^e) (icfe) I 1 "?)- (7) 

The dimensionless constants ci,2 are quoted in Table 1 for fi r = 0.5, 0.1, 10~ 2 , 10~ 3 , and 10~ 4 . For all /i r , the star is oblate 
for A > 0.4 and prolate for A < 0.4. The functional dependence of e on A stays roughly the same as p r changes. The smaller 
magnetic permeability of the superconducting stellar matter enhances the density perturbation by a factor of ~ fj,^ , as 
evident from the force balance equation [Eq. (3)], and this is embodied in the approximate scaling ci oc /i^T 1 for small y, r . 

We can now ask how much closer the theoretical curves approach the data when superconductivity is included. Looking 
at Fig. [2] for example, the curve for (say) p r = 10 -3 and A = 10 -3 (not drawn) is higher than the lower black dashed curve 
by three orders of magnitude, which is still below observational limits. We conclude that superconducting interiors are easily 
compatible with current observational upper limits, if the external magnetic field is not reduced by accretion in any way. 

We caution that, in this first pass, we assume fj, r ^ 1 throughout the star, instead of only in some region. This assumption 
is physically implausible and is only taken to simplify our (illustrative) calculation. A more thorough calculation where fi T is 
allowed to vary inside the star is needed before definite conclusions can be drawn regarding the effects of core superconductivity 
on ellipticity. However, we conjecture that the calculation presented in this section sets the upper limit on the effects of a 
superconducting interior on stellar deformation: if superconductivity is limited to a smaller region in the star, the changes to 
e (relative to fi = fig) will be less than predicted by Eq. (7) and Table 1. 



3.2 Accretion 

Recycled pulsars are selected as the subjects of this study for two reasons: (i) their spin-down rates are lower than those of 
other objects, yielding more stringent indirect gravitational-wave limits; and (ii) their internal magnetic fields may be much 
greater than their surface fields due to diamagnetic screening or burial, yielding larger hydromagnetic deformations than one 
might otherwise expect. We now examine point (ii). Let us ask what happens if the actual magnetic field strength just below 
the surface takes its pre-accretion value (e.g., before diamagnetic screening or burial) -B s , a ctuai = -B s , observed /(,, where £ is some 
dimensionless 'shielding factor'. We note in passing that this scenario is more realistic than those considered in Sec. 3.1: it is 
unlikely that the core poloidal fiel d is reduced to ~ 10 4 T in a rec ycled pulsar like the surface field, given the high electrical 
conductivity expected in the core l|Goldreich fc Reisenegger 1992t ). except in the special situation where the source currents 



reside exclusively in the crust. 

We recalculate e using Eq. (7), substituting B s , actual for B B and keeping fi r = 10 -3 . We plot the results for £ = 10 -4 as 
the red curves for A = 1 (bottom solid red curve) and A = 10 -3 (top solid red curve) in Fig. [2] We also plot e for the case 
of £ = 10~ 2 as the blue curves for A = 1 (bottom solid blue curve) and A = 10~ 3 (top solid blue curve). For clarity, the 
case fj, — fxo is not presented; it lies three decades lower than the curves with /i r = 10 -3 . Without a better knowledge of the 
screening/burial process than is c urrently at hand, it is best simply to bracket the plausible range 10~ 4 < £ ^ 10" 2 inferred 
from population synthesis studies I Kiel et al. 20081 ). 



Using the potentially stronger shielded pre-accretion fields, the e(-B s ) curves in Fig. [2] come close to the spin-down limits 
on e. Now we can see that, aside from perhaps one object, namely PSR J 1823-3021A, all the pulsars plotted must have some 
internal toroidal field component, i.e. A < 1, if one has £ ^ 10 -2 . From the red curves and the dots in Fig. [2l many objects 
seem to have A < 0.01, even A < 10 -3 in the notable case of PSR J1910-5959C. On the other hand, the heavily-shielded case 
A = 10" 3 , ^ = 10~ 4 (for example) seems to be ruled out by observations (top solid red curve in Fig. [2]). In general, A cannot 
be too small, for several reasons: it is ruled out by observations, the virial limit sets an absolute upper bound on the internal 
field strength, and small A leads to an unstable field configuration (|Braithwaitell2009[). By contrast, large A is not ruled out 
by observations; the only upper bound (A = 0.8) is set by the stability analysis of iBraithwaitd |2009j). Fig. [2] also tells us that 
^ cannot be too small, otherwise the solid red curves exceed the observational upper limits (for £ sC 10~ 4 ). 
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(b) 



0.0 0.5 1.0 1.5 2.0 




N 0.2 - 





Figure 3. Poloidal field lines inside and outside a superconducting star with fi r = 10 _1 [(a) and (b)] and fir = 10 — 3 [(c) and (d)]. The 
surface of the star is represented by the dashed hemisphere. The right panes (b) and (d) zoom in to the boxed region around the neutral 
circle in the left panels (a) and (c) respectively. Note that all field lines are, in fact, closed, since the magnetic field is divergence- free. 



The data in Fig. [2j together with the stability-based limits on A set by Braithwaitel ( 20091 ) . allow us to use gravitational 
wave observations and measurements of surface fields to infer bounds on both A and the shielded field -B s , actual, in principle. 
We show several possibilities in Fig. 0] where we draw curves of constant e = 1.184 x 10~ 9 (corresponding to the current 
lowest spin-down upper limit, see Fig. [2| for a given value of -B s , observed (solid: 10 5 T, dashed: 10 4 T), for fi r — 1 (blue curves) 
and fi r = 10~ 3 (red curves). To be consistent with spin-down measurement, an object must lie above the curve relevant to 
the scenario being considered. The shaded region bordered by thick black dashed lines indicates the t heoretical limits on A 
[from stability; 10" 3 < A < 0.8 l|Braithwaitell2009h ] and £ [from population synthesis; 10" 4 < £ < 10" 2 jKiel et alj|200s| )]. In 



other words, Fig. [?] allows us to determine the allowed combinations of A and £ for a given observed surface field and e limit 
(here e = 1.184 x 10 -9 ). For example, a detection from a pulsar with surface field of 10 4 T indicates 10 -3 ^ A ^ 7 x 10~ 3 
and 10~ 4 < £ < 2.5 x 10~ 4 , if we assume the star is non-superconducting (the blue dashed curve), or 10" 3 < A sC 4 x 10 _1 
and 10" 4 ^^8x 10 -3 (prolate star) or 4.5 x 10" 1 sC A sC 0.8 and 10~ 4 < £ sC 2.4 x 10" 4 (oblate star) if the star is 
superconducting with fi r = 10 -3 (the red dashed curve). From Eq. (0, we see that the oblate cases (with ci ^ A ^ 1) are 
weakly deformed relative to the prolate cases (10 -3 ^ A ^ c^). Hence, for oblate stars to be readily detectable, we need small 
£ and /i r to boost e (red curves). 



4 DISCUSSION 



This paper makes a first attempt at combining non-barotropic magnetized stellar models (Mas trano et alJl201lT ) with LIGO 
non-detections and radio timing data to constrain the ra ti o A and hence the inte rnal magnetic field of recycled pu l sars. The 
models resemble the 'twisted torus' of Ciolfi et al. ( 2009h . Lander h, Jones! ( 20091 ') . and Ciolfi. Ferrari, fc Gualtieri 1 2010l ). in 
that they consist of a potentially strong internal toroidal field, which is not observable directly and is confined inside the 
star, and an ex ternal field, which can be measured. The main difference is that we assume the star is stably stratified but 
non-barotropic (Reiseneggcr 2009 1 ). Because of this, our models allow the poloidal and toroidal components to be adjusted 
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Figure 4. Curves of constant ellipticity (t = 1.184 X 10 -9 ; strictest limit from Fig. [2} on the A— £ plane, where £ is a dimensionless 
'shielding factor', £ = B s bserved/-Bs, actual) an d A is the ratio of the internal poloidal magnetic field energy to the total internal magnetic 
field energy. Blue curves: non-superconducting fi r = 1, with B s .observed — 

10 4 T (dashed) and B s , 

observed — 10^ T (solid). Red curves: 

superconducting /j, r = 10~ 3 , with -B s , observed = 10 4 T (dashe d) and B s , observed = 10 5 T (solid). The thick black dashed li nes are the 
theor etical bounds on A and £: from stability 10 -3 ^ A ^ 0.8 (Braithwaite 200^) and population synthesis 10 -4 ^ £ ^ 10' " 2 (iKiel etajj 
1200^1 . 



independently of each other; in barotropic models (such as the twisted torus model), the poloid al and toroidal field s must 
obey certain relations to ensure that pressure is always a function of density purely [see also, e.g.. lHaskell et al.1 (j2008l >] . This 



means that our models c an easily accommodate fi e ld con figurations with 0.01 < A < 0.8, the range found to be stable by the 
numerical simulations of Braithwaite fc Nordlundl ( 20061 ) . 

As seen in Fig. [2] our calculated e is much lower than the LIGO and spin-down upper limits when the inferred dipole 
field strengths are input directly into Eq. (5)0 Note that the lower black dashed curve in Fig. [2] already assumes a strong 
toroidal field (A = 10~ 3 ); a weaker toroidal field will generate a curve even lower down. Next, we calculate how interior 
superconductivity (jj, r — 10 -3 ) and accretion-induced screening or burial (red and blue bands in Fig. [2]) enhance e. This latter 
possibility is further explored in Fig. 0] which shows how, in principle, a gravitational wave detection from a recycled pulsar 
with a certain surface field (measured from spin down) can be used to infer both A and shielding factor £. Fig. [4] also tells us 
the sets of parameters that are ruled out by current observations: for a given pulsar with a measured spin-down field strength, 
non-detection means that £ and A lie above that particular curve in Fig. [4] For example, the upper solid curve of Fig. [2] 
corresponding to £ = 10~ 4 and A = 10 -3 , lies at the corner of the theoretically allowed shaded region in Fig.|4]and is therefore 
ruled out even for a non-superconducting star with a relatively weak observed field strength of 10 4 T (because it lies below 
the dashed blue curve in Fig. |4]). 

I s it pos s ible th at t he magnetic fi eld of a recycled pulsar is much stronger just below the surface than the inferred dipole 
field? iKrolikl |l99ll ) and UronsI jl993l ). motivated by the discrepancy between field strengths inferred from dipole spin down 
and from cyclotron line measurements of a ccretin g X-ra y pulsars, raised such a possibility. They proposed that higher-order 
multipoles may exist close to the surface. Aronsl (1993) then showed that surface field strengths ~ 10 4 times higher than 
the observed dipole field can account for the aforementioned discrepancy, as well as the location of the millisecond pulsars 
on the P-P diagram. The surface field can also be masked by magnetic field burial, whereby accreted matter compresses 
the polar magnetic flux into a narrow be lt around the equator. It has been shown that accretion of > 10~ 5 Mp) is enough 
to alter the dipole mom ent significantly (Lai 19991 : Payne fc Melatos 2004 : Melatos fc Pavnei 20051 : Zhang fc Koiima 2006; 



Vigelius fc MelatosH2009T). d epending on the equation of state IjPrivmak. Melatos. fc Pavnei 1 201 if ). 



Pons fc Gepperti (|2007l ) raised the possibility that the large-scale magnetic field of a neutron star is supported by long- 



lived currents in its superconducting core, while small-scale, fast-decaying (decay time 10 5 -10 7 yr) magnetic structures also 



4 Recall that the magnetic axis is assumed to be perpendicular to the rotation axis throughout this paper. 
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exist near the s urface, supported by sho rt-lived currents in the crust. After trying several different initial poloidal and toroidal 
field strengths, iPons fc Geppertl (2007) found that all their models eventually (after ~ f Myr) reconfigure into a long-term 
stable configuration comprising a dipolar poloidal component and a quadrupolar/octupolar toroidal component. The decay 
of the crustal currents is domin a ted by Hall drift in the first 10 3 -10 4 yr if the initial fi eld is > 10 10 T, i.e. magnetar strength 
( Goldreich fc Reisenegger] Il992l ; Rheinhardt fc Geppertl I2OO2I ; Pons fc Geppertl |2007| ) . When the external dipole field is of 
ma gnetar strength, the inte rnal crustal field can still be ~ 10 times greater. 

Pons fc Geppert 12007) did not calculate the magnetic deformation of their star, but one can extend their work and 



that of iMastrano et all ( 2011 ') to investigate the effects of this strong, internal, quadrupolar/octupolar toroidal field on 
stellar ellipticity. To do so, one must remember that Eq. (5) applies to an internal field that is continuous with an external 
di pole, as shown in Fig . IT] W hen t he field is distorted signi ficantly by accretion [into, e.g., the 'equatorial tutu' shape found 
by iMelatos fc Phinnevf (|200ll ) and IPavne fc Melatosl (|2004h ]. then the internal toroidal field described by Eq. (2) must be 
recalculated to ensure that the total field still has the desired properties listed at the start of Sec. 2. A similar adjustment 
must be carried out to accommodate multipolar external fields. It can be shown, for example, that one needs t wo toroidal 
field belts around the new neutral circles at (0.7R, n/4) and (0.7R, 37r/4) to match a quadrupolar external field (Mastrano 
2010h , 

Lastly, we remind the reader that a strong internal field is not the only possible cause for stellar deformation. Other 
effects that can deform a star are: 

• rotat ion; ellipticity induced by centrifuga l forces is en ~ 0.3(V/kHz) 2 , where v is the spin frequency of the neutron 
star jPines fc Shahamlll972l ; ICutlerl 12002 ). Magnetar spin frequencies are ~ 0.1 Hz, giving eq ~ 10 9 . Ho wever, the 
defor mation is aligned with the rotation axis [except for a fraction ~ 10~ 4 prop ortional to the shear modulus ( Goldreich 
Il970l )]. so it does not contribute to gravitational wave emission or precession jMelatoslboOOri 



^sTal shear stresses; ellipticity supported by the crust's <-I;im icii v ill biskoll. .Tones, fc Anderssonl B. up to a maxi 



mum of ~ 10~ for conventio nal neutron stars, ~ 10~ for hybrid q uark-baryon or meson-condensate stars, or ~ 10~ 
for solid strange quark stars (|Owe n 2005; H orowitz fc Kad au 2009). 

nuclear processes; a sustained, asymmetric temperat ure step of ~ 5% at the base of the cr ust drives electron capture 
reactions which produce compositional asymmetries ( Ushomirskv. Cutler, fc Bildsten 2000h . 
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